Response properties of III-V dilute magnetic semiconductors: interplay of disorder, 
dynamical electron-electron interactions and band-structure effects 
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A theory of the electronic response in spin and charge disordered media is developed with the 
particular aim to describe III-V dilute magnetic semiconductors like Gai-^^Mua^As. The theory com- 
bines a detailed k ■ p description of the valence band, in which the itinerant carriers are assumed 
to reside, with first-principles calculations of disorder contributions using an equation-of-motion ap- 
proach for the current response function. A fully dynamic treatment of electron-electron interaction 
is achieved by means of time-dependent density functional theory. It is found that collective excita- 
tions within the valence band significantly increase the carrier relaxation rate by providing effective 
channels for momentum relaxation. This modification of the relaxation rate, however, only has a 
minor impact on the infrared optical conductivity in Gai-^^Mua; As, which is mostly determined by 
the details of the valence band structure and found to be in agreement with experiment. 

PACS numbers; 72.80Ey, 75.50Pp, 78.20.Bh 
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I. INTRODUCTION 

The idea of using both charge and spin of electrons 
in a new generation of electronic devices constitutes the 
basis of spintronics.>i The magnetic properties of the ma- 
terial combined with its semiconducting nature makes 
dilute magnetic semiconductors (DMSs) potentially ap- 
pealing for various spintronics applications. - In particu- 
lar, the effect of carrier mediated ferromagnetism opens 
up the possibility to control the electron spin and mag- 
netic state of a system or device by means of an electric 
field. A lot of attention is drawn to Gai-j^Mn^jAs due to 
the well developed technology of the conventional GaAs 
based electronics and discovery of its relatively high ferro- 
magnetic transition temperaturCj— with a current record 
of Te = 185 K,^ 

Unlike most other III-V DMSs, the nature of the itin- 
erant carriers in Gai_a:Mna;As is still under debate4ii^ It 
is widely accepted that for low-doped insulating samples 
the Fermi energy lies in a narrow impurity band. For 
more heavily doped, high-Tj, metallic samples there are 
strong indications that the impurity band merges with 
the host semiconductor valence band forming mostly 
host-like states at the Fermi energy with some low-energy 
tail of disorder-related localized statesi^ First-principles 
calculations2r^ have so far not been fully conclusive re- 
garding the nature of the itinerant carriers in this regime, 
and further theoretical studies continue to be necessary. 
Meanwhile, attention has shifted to model Hamiltonian 
approaches assuming either the valence band^*^ or impu- 
rity band^^ picture and their ability to adequately de- 
scribe the experimental results in Gai_a;Mna;As. 

Most theoretical approaches assuming the valence- 
band nature of itinerant holes in Gai-^^Mn^^As treat the 
band structure in detail, while disorder and many-body 
effects are only accounted for using simple phenomeno- 
logical relaxation time approximations and static screen- 
ing modelsJ^^— On the other hand, the extreme sensitiv- 
ity of magnetic and transport properties of Gai_2;Mn2:As 



to details of the growth conditions'^ and post-growth 
annealingi^"— points to the crucial role played by the 
defects and their configurations, and has stimulated in- 
tense research on the structure of defects and their in- 
fiuence on the various properties of the systemji^ It is 
essential, therefore, to develop a theory of electrical con- 
ductivity in DMSs with more emphasis given to disorder 
and electron-electron interactions, without neglecting the 
intricacies of the electronic band structure. 

Here we present a comprehensive theory for the elec- 
tron dynamics in DMSs which accounts for the complex- 
ity of the valence band structure of the semiconductor 
host material and treats disorder and electron-electron 
interaction on an equal footing. In previous work we 
used a simplified treatment of the semiconductor va- 
lence band^"'^^ or considered only static properties of the 
system. In this paper we simultaneously account for the 
complexity of the valence band, use a first-principles ap- 
proach to describe disorder contributions, and employ a 
fully dynamic treatment of electron interactions. 

To account for the valence band structure we use the 
generalized k • p approaches where a certain number of 
bands are treated exactly while the contribution from the 
remote bands is included up to second order in momen- 
tum. To describe disorder effects we use the equation of 
motion for the paramagnetic current response function 
of the fully disordered system. This approach has some 
similarities to models developed earlier using the memory 
function formalism. ^4"™ The advantage of our approach 
as compared to the memory function formalism is the rel- 
ative simplicity and transparency of the derivation and 
the straightforward possibility to include the spin degree 
of freedom. Another advantage is that our formalism is 
expressed in terms of a current-current and a set of den- 
sity and spin-density response functions. This enables us 
to use the powerful apparatus of time-dependent density- 
functional theory (TDDFT)^^ to treat many-body effects 
such as dynamic screening and collective excitations of 
the itinerant carriers in principle exactly. 
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The paper is divided into two major sections and con- 
clusions. For ease of reading, some of the derivations are 
presented in appendices. The theory section (Sec. HIl) 
is organized as follows. In Sec. Ill Al we present our 
general formalism based on the equation of motion of 
the current-current response function of the disordered 
system. In Sec. Ill Bl we describe the evaluation of the 
current-current, density and spin-density response func- 
tions for the multiband system using a generalized k • p 
perturbation approach. Next, in Sec. Ill CI we show the 
treatment of electron-electron interaction by means of 
TDDFT. In Section IIIII we first discuss the new fea- 
tures that the valence band character of itinerant car- 
riers brings into the system, namely the dominance of 
the long- wavelength side of the single-particle excitation 
spectrum by the interband spin transitions and the ef- 
fective suppression of the collective plasmon excitations 
within the valence band for the whole range of momen- 
tum. Next, in Sec. IIII Bl we discuss the effect of mag- 
netic doping: spin and charge disorder in the system and 
modification of the band structure in the magnetically 
ordered phase. We show that the full dynamic treatment 
of electron-electron interactions allows us to capture the 
effect of collective excitations on the carrier relaxation 
time. We then compare our results also with experimen- 
tal data on infrared conductivity. Finally, in Sec. IIVI we 
draw our conclusions. 



II. THEORY 
A. General formalism 

We discuss a system described by the Hamiltonian 

H = He+ H„r + Hd, (1) 

where He is the contribution of the itinerant carriers and 
Hm represents the subsystem of localized magnetic spins. 
These two terms constitute the "clean" part of the total 
Hamiltonian. The last term in Eq. ([ij describes disorder 
in the system: 



with the components 
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where the four-component charge and spin disorder scat- 
tering potential 
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is coupled to the four-component vector of charge and 
spin density operators of the itinerant carriers: 
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Here, cr^ (/i = ) is defined via the Pauli matri- 

ces, where is the 2x2 unit matrix, cr* = (cr^ ±i(T^)/2, 
and |wn,q) are the two-component Bloch function spinors 
with wave vector q and band index n. The summation in 
Eq. (jS)) is performed over all defects. Note that the mean 
field part of the p-d exchange interaction between itiner- 
ant holes and localized spins is absorbed into the clean 
system band structure Hamiltonian iJg; disorder in our 
model consists of the Coulomb potential of charge defects 
and fiuctuations of localized spins around the mean field 
value {S). 

The general case of multiple types of defects, includ- 
ing defect correlations, was considered in Ref. [20. For 
simplicity we here include only the most important de- 
fect type, namely randomly distributed manganese ions 
in gallium substitutional positions (Mnca)- Our model 
treats localized spins as quantum mechanical operators 
coupled to the band carriers via a contact Heisenberg in- 
teraction featuring a momentum-independent exchange 
constant J. We use the value of V J — —55 meVnm"^, 
which corresponds to the widely used DMS p-d exchange 
constant NqP = — 1.2eVJ^ The z-axis is chosen along 
the direction of the macroscopic magnetization. 

Earlier we developed a theory of transport in charge 
and spin disordered media with emphasis on a treatment 
of disorder and electron-electron interaction<2i It is based 
on an equation of motio n^^i^^ approach for the param- 
agnetic current-current response of the full, disordered 
system: 

z 

Xjp„ip^(r,r',T) = --e{T){[jpo,{T,r),jpi3{r')])H, (6) 



where 



(7) 



is the paramagnetic current-density operator in Heisen- 
berg representation and a, P = x,y, z are Cartesian coor- 
dinates. During the derivation we assumed our system to 
be macroscopically homogeneous, which is justified if the 
coherence length of the electrons is much shorter than 
the system size. In this case, summing over all electrons 
will leave us with an averaged effect of disorder that does 
not depend on the particular disorder configuration. For 
such macroscopically homogeneous systems the response 
at point r depends only on the distance |r — r'| to the per- 
turbation and not on the particular choice of points r and 
r'. Another major approximation involves the decoupling 
procedure, where we neglect the influence of the itinerant 
carriers on the localized spins. Therefore, our approach 
does not include magnetic polaron effects and lacks the 
microscopic features of carrier mediated ferromagnetism. 
The latter, however, can be reinstated to some extent by 
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introducing a phenomenological Heisenberg-like term in 
the magnetic subsystem Hamiltonian . Details of the 
derivation are presented in Ref. lU 

The final expression for the total current response 
reads 
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where Xpi^p'^O^T^) is the set of charge and spin density 
response functions with respect to the operators (|4])-(l5]) 
and the superscript "c" indicates quantities defined in 
the clean system. By comparing Eq. ([8]) with the Drude 
formula in the weak disorder limit cut ^ 1, 
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we identify the tensor of Drude-like frequency- and 
momentum-dependent relaxation rates of the form 



k 

X (xp-p" (q - k, c^) - Xp.p. (k, 0)) . (10) 



Note that the right-hand side of Eqs. (|8]) and (|T0|) con- 
tains the set of spin and charge response functions of 
the full, disordered system. Therefore, strictly speak- 
ing, Eq. ((HI) should be evaluated self-consistently^° with 
the continuity equations closing the loop. Here we use a 
simplified approach based on two approximations. First, 
taking the weak disorder limit in the right hand side of 
Eq. (fTO]) we retain terms up to the second order in com- 
ponents of the disorder potential. In other words, the 
spin and charge response functions of the full system in 
Eq. (jlOl) are replaced by their clean system counterparts: 



XpA'p-(q- k,w) 



(11) 



Next we assume that the paramagnetic current response 
function of the full system may be expressed as the 
clean system response function with a lifetime broaden- 
ing given by Eq. (fTO| : 



(q' ^) ~ (q, ^ - (12) 



Equations pO)) ~(IT2 |) will be used in the following. 



B. Multiband k ■ p approach 

In order to obtain the conductivity through Eqs. (fTO|)- 
([T2l) we will have to calculate the paramagnetic current 
response and spin and charge density response functions 
of the clean system. To properly describe the complex- 
ity of the semiconductor valence band we are going to 
implement the multiband k • p approach. 

First we derive the current and density response func- 
tions in the formal basis of the Bloch states 
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Within second quantization in the basis (jl3p , the param- 
agnetic current in the system with spin-orbit interaction 
is given by 



which diagonalize the clean system Hamiltonian 
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with 



TT = p 



[<T X VC/e], 



(16) 



where Uc is the periodic crystal field potential. Hereafter, 
performing the real space integration we assume that the 
envelop function varies slowly on the scale of the unit 
cell. 

Introducing the time dependence of the creation and 
destruction operators in (fTSj) . the paramagnetic current 
response of the multiband system can be directly evalu- 
ated, and one finds 
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A similar procedure for the spin and charge density response yields: 
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All we need now for evaluating Eqs. ([T7)) and (|18p is to determine the form of the periodic Bloch functions |u„^k) that 
diagonalize the clean system Hamiltonian. The common approach is to diagonalize the multiband k • p Hamiltonian 
that treats certain bands exactly and treats contributions from remote bands up to second order in momentum. The 
derivation of such a Hamiltonian is outlined in AppendixlX] By diagonalizing the matrix of this Hamiltonian, however, 
we obtain the eigenvectors of the modified Hamiltonian (|A7[) . Before evaluating the matrix elements between Bloch 
periodic functions |wn_k) in Eqs. p7)) and (|18l) we therefore have to perform the unitary transformation (jA4[) . Details 
of these calculations are presented in Appendix |B1 

The final expression for the paramagnetic current response function in the long-wave limit q = (since we are 
looking for the optical response) is given by 



^ n,n ,k 
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where H denotes the effective multiband k- p Hamiltonian (IA7|) and B(n, k) is its eigenvector for the state with 
energy e„^k- The charge and spin density response is approximated by 
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If o-^ = (5-'^)+, i.e. for Xnn, Xs^s^ and Xs±s=Fj the second sum is a real quantity. Then, the imaginary part is 



^^[Xpf.(pf.)+(q,w)] = - J d^Hfn'M-q - /«,k)5[?iw - {en,k " £n',k-q)] 
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It is seen that in the long-wavelength limit (q — > 0) the 
imaginary part of the density response [p^ = tr^) van- 
ishes as a product of orthogonal states, while the imag- 
inary part of spin response is, in general, finite. We 
conclude from this that the long-wavelength spectrum 
of single-particle excitations is dominated by spin transi- 
tions. 

The calculations were performed within an 8-band 
k • p model. The basis functions and explicit form of 
the Hamiltonian matrix are presented in Appendix [Cl 



C. Electron-electron interaction 

A major advantage of our formalism is that it is ex- 
pressed in terms of current and density response func- 
tions. This allows us to use the powerful apparatus of 
TDDFT to account for the effects of electron-electron 
interaction. 

Let us first examine the current response of the clean 



system. In this paper we are considering the optical re- 
sponse, i.e. the response to transverse perturbations. 
Since transverse perturbations only induce a transverse 
response in a homogeneous system, there are no density 
fiuctuations directly created by an electromagnetic field. 
The total current response of the interacting system in 
this case can be expressed as 

( 1 f J 47re cj^ 

(x (q,^)j =(xo(q>w)j + ^2^g2 2 + — ^9Gt+, 

(22) 

where xi is the response of the noninteracting system, Vq 
is the Coulomb interaction, and the local field factor Gt+ 
represents corrections from the exchange-correlation (xc) 
part of the electron interaction. 

The corrections to the transverse current response 
function caused by electron-electron interaction are rela- 
tivistically small in this case and can be neglected. So, 
for the transverse current response of the clean system 
we will use the noninteracting form. 

The set of the density and spin-density response func- 
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tions of the clean system enters our expression (jTOj) 
for the frequency- and momentum-dependent relaxation 
rates. TDDFT allows us to describe all the effect of 
electron interaction, including correlations and collective 
modes, in principle, exactly. Within the TDDFT formal- 
ism the charge- and spin-density responses of the inter- 
acting system can be expressed as)^ 



Xo 



-^q,uj)-u{q)-Lci^,^), (23) 



where all quantities are 4x4 matrices and xo denotes 

the matrix of response functions of the noninteracting 
system, u{q) is the Hartree part of the electron-electron 

interactions, and /xc represents xc corrections in the form 

of local field factors. As a simplification we use only the 
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exchange part of /xc and apply the adiabatic local spin 

density approximation. Explicit expressions for the local 
field factors of the partially spin polarized system are 
given in Appendix IdI 

In general, /xc is a symmetric 4x4 matrix. If, however, 

the z-axis is directed along the average spin, then the 
ground state transversal spin densities vanish, — p- — 
0, and the matrix /xc becomes block-diagonal: 
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Performing the matrix inversion in Eq. ()23|) we obtain the 
tensor of response functions of the interacting system in 
the form 
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ELFF = 1 - {v{q) + /ll)xL(9,w) - fzzX%s4q,^) - flz (Xns^q,^) + Xs^nil,^)) + {fzz{v{q) + /ll) - fl) A, (26) 



and 



A. Clean p-type GaAs 
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III. RESULTS AND DISCUSSION 

We now discuss applications of our formalism for the 
specific case of GaMnAs DMSs. The band structure pa- 
rameters used in our calculations correspond to those of 
the GaAs host material: the band gap and spin-orbit 
splitting are Eg = 1.519 eV and A = 0.341 eV, Luttinger 
parameters are 71 = 6.97, 72 = 2.25 and 73 = 2.85, con- 
duction band effective mass is nie — 0.065 mo, Kane mo- 
mentum matrix element is Ep — 27.86 eV and the static 
dielectric constant K = 13. The s{p)-d exchange interac- 
tion constants within the conduction and valence bands 
are N^a = 0.2 eV and Nq[3 = —1.2 eV, respectively. 



Before considering the effects of magnetic impurities 
and associated charge and spin disorder on the transport 
properties, we would like to discuss some new features 
that the valence band character of the itinerant carriers 
brings into the system. They stem from the complex- 
ity of the semiconductor valence band: strong spin orbit 
interaction and the F-point degeneracy of the p-states. 

The multiband nature of the valence band gives rise 
to a rich single-particle excitation spectrum. In Fig. [1] 
we show a schematic representation of the valence band 
structure of a p-type semiconductor. Arrows indicate the 
possible single-particle excitations. In addition to the in- 
traband excitations within the heavy hole band (anal- 
ogous to the excitations within the conduction band of 
n-doped semiconductors), here we have intra-band exci- 
tations within the light hole band as well as inter- valence 
band excitations between light and heavy hole bands and 
between split-off and heavy and/or light hole bands. 



FIG. 1. (Color online) Schematic diagram of the possible 
single-particle excitations in the valence band of a p-type 
semiconductor. Dashed lines indicate intra- valence band exci- 
tations within the heavy hole band (hh), within the light hole 
band (Ih) and inter-valence band excitations between heave 
hole and light hole bands (hh-lh) and between split-off and 
heavy hole and light hole bands (so). 
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FIG. 2. (Color online) Imaginary part of the noninteracting 
density and longitudinal spin response functions in p-doped 
GaAs for different wave vectors q = O.OOSA (a) and q = 
O.OSA"^ (b). The hole concentration is p = 3.5 x 10^° cm^^. 

The variety of the possible single-particle excitations 
substantially modifies the density and spin response of 
the system. Some of the modifications are not very obvi- 
ous. At the end of the Sec. IIIBI we already mentioned the 
significant difference between spin and density responses 
in the long wavelength limit. Let us consider this in more 
detail. The spin response of the noninteracting electron 
gas coincides with the density response and can be ex- 
pressed through the Lindhard function. The spin-orbit 
interaction within the valence band breaks down this cor- 
respondence. 

In Fig. [5] we plot the imaginary part of the noninter- 
acting density and longitudinal spin response functions 
in p-doped GaAs for different wave vectors. For a small 



FIG. 3. (Color online) Schematic diagram of the excitation 
spectrum within the semiconductor valence band. Labels in- 
dicate the edges of single-particle excitation regions within 
the heavy hole band (hh), within the light hole band (Ih), 
between heavy hole and light hole bands (hh-lh) and between 
the split-off band and heavy and light hole bands (so), see 
Fig. [T] As a result, the plasmon mode in the valence band 
lies entirely within the single-particle excitation spectrum and 
is effectively suppressed due to Landau damping. 



wave vector q — O.OOSA the longitudinal spin response 
exhibits a strong peak around 0.2 eV associated with 
inter-valence band spin excitations between heavy and 
light hole subbands. The corresponding density excita- 
tions are suppressed due to the orthogonality of the initial 
and final states, see Eq. (I?T|) . As a result, the density re- 
sponse for short wavevectors is almost nonexistent. If we 

increase the wave vector to q — O.OSA , the intraband 
excitations within the heavy hole band become notice- 
able in both density and spin responses. The longitudi- 
nal spin response, however, still prevails in the range of 
inter-valence band transitions. 

This leads us to conclude that the long-wavelength 
spectrum of the single-particle excitations in p-doped 
semiconductors is dominated by the inter-valence band 
spin excitations. The origin of this effect is in the spin- 
orbit interaction, which mixes spin and orbital degrees 
of freedom. Without the spin-orbit interaction, vertical 
spin excitations would be prohibited due to the orthogo- 
nality of the orbital parts of Bloch functions. 

Another interesting feature of p-doped semiconductors 
is the effective suppression of the collective modes in the 
valence band. In the conventional picture of the con- 
duction band, collective plasmon excitations are well de- 
fined in the long-wavelength side of the excitation spec- 
trum. With increasing momentum, the collective mode 
approaches and then enters the region of single-particle 
excitations, where it becomes rapidly suppressed due to 
Landau damping. 




FIG. 4. (Color online) Dispersion (dashed black) and lifetime 
broadening (solid red) of the valence band plasmon calculated 
for the p-doped GaAs with the hole concentration of p — 
3.5 X 10^" cm~'^. Dotted lines correspond to the onset of the 
intraband single-particle excitations within the light hole and 
heavy hole bands. 



The situation is different for the valence band. In 
Fig. [3] we plot a schematic diagram of the excitation spec- 
trum. The excitation region for single-particle transitions 
within the heavy hole band is qualitatively similar to that 
of the conduction band. In the valence band, however, 
the single-particle excitation spectrum is extended due to 
the intraband transitions within the light hole band and 
interband transitions between heavy and light hole bands 
and between split-off and heavy /light hole bands (red 
and blue arrows in Fig. [!}. In Fig. [3] the corresponding 
regions of single-particle excitations are shaded with dif- 
ferent patterns. It can be seen that the collective mode in 
the valence band falls entirely within the region of single- 
particle excitations and, therefore, becomes suppressed 
even at the long-wavelength side of the spectrum. Error 
bars in Fig. [3] indicate the plasmon resonance broadening 
due to Landau damping. 

To illustrate the effect we have performed numerical 
calculations of the plasmon dispersion and the lifetime 
broadening of the collective excitations in the valence 
band of the p-doped GaAs. The plasmon frequencies 
were determined as the zeros of the real part of the RPA 
dielectric function and the lifetime broadening is asso- 
ciated with the imaginary part of the frequency poles. 
In Fig. U the black and red lines correspond to the dis- 
persion and the lifetime of the plasmon excitations, re- 
spectively. The dotted lines indicate the regions of the 
intraband single-particle excitations within the light hole 
and heavy hole bands, compare with Fig. [31 At small 
wavevectors the plasmon mode falls within the region of 
inter-valence band single-particle excitations resulting in 
a lifetime broadening of the collective resonance of about 
5 meV. Once the plasmon dispersion enters the region of 
single-particle excitations within the light hole band, the 
life-time broadening substantially increases into the 30- 



FIG. 5. (Color online) Total (charge and spin) carrier re- 
laxation rate for Gao.948Mno.o52As with hole concentration 
p = 3 X 10^" cm"'^. Dashed line: static screening model. 
Solid line: evaluation of Eq. ^ with full dynamic TDDFT 
treatment of electron interaction. See discussion in text. 
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FIG. 6. (Color online) Infrared conductivity of ferromagnetic 
Gao.g48Mno.o52 As with hole concentration p = 3 x 10^'^ cm~^. 
Calculations are performed according to Eq. (|12|l . and using 
a relaxation rate obtained through Eq. (|lUp (solid line) or a 
fixed = 230 ps~^ (dashed line). Symbols are the experi- 
mental data of Ref. [SJ. 

40 meV range. An additional sharp rise in the damping 
takes place when the collective mode enters the region of 
heavy hole intraband excitations. 



B. Magnetically doped GaMnAs 

The introduction of magnetic impurities in GaAs has 
two consequences. First, charge and spin disorder are 
brough into the system and, second, the mean field part 
of the p-d exchange interaction between localized spins 
and itinerant holes causes modifications of the valence 
band structure once the system enters the magnetically 
ordered phase. 
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FIG. 7. (Color online) Band structure of ferromagnetic 
Gao.95Mno.05As with hole concentration p = 3.5 x 10^° cm~^. 
Alignment of localized spins results in strongly anisotropic 
valence band spin splitting. Inset shows a cut of the Fermi 
surface by the plane ky = 0. 



Let us consider the effect of disorder first. In calcu- 
lating carrier relaxation rates, most theoretical models 
for GaMnAs use a static screening approach, where all 
many-body effects are reduced to the static screening of 
the Coulomb disorder potential. Within our model, how- 
ever, the momentum and frequency dependent relaxation 
rate pop is expressed through the set of density and spin- 
density response functions that allows us to use the full 
dynamic treatment of electron-electron interaction, thus 
accounting for the variety of many-body effects including 
correlations and collective modes. 

In Fig. [5] we plot the frequency dependence of the 
total (charge and spin) relaxation rate calculated for 
Gao.948Mno.o52As within the static screening model and 
using the full dynamic treatment of electron-electron in- 
teraction according to Eq. ([TU]). The difference between 
the two curves in the static limit is due to the xc part of 
the electron-electron interaction that affects both charge 
and spin scattering. The most striking difference, how- 
ever, is the pronounced feature appearing between 0.2 
eV and 0.5 eV associated with the collective modes. Al- 
though we have seen above that the collective excitations 
are significantly damped in the valence band, they still 
play an important role in the transport properties of the 
system providing an effective channel for momentum re- 
laxation. Their contributions give up to 50% increase 
to the total carrier relaxation rate. Note that, due to 
their longitudinal character, the plasmon modes do not 
directly affect the optical response and enter only indi- 
rectly through the tensor of frequency and momentum 
dependent relaxation rates ([TU]). 

In Fig. [B] we compare our calculations of the infrared 
conductivity of ferromagnetic Gao.948Mno.052 As with the 
experimental data of Singley et ali^ The calculations 




u [eV] 

FIG. 8. (Color online) Temperature dependence of infrared 
conductivity of Gao.948Mno.052 As with hole concentration p = 
3 X 10^° cm~'^ and Tc — 70 K calculated with weak disorder, 
lifetime broadening of F = 5 meV. 

were performed according to Eq. p^ . Solid line corre- 
sponds to a relaxation rate obtained through Eq. (|10l) . 
dashed line describes calculations with the fixed = 
230 ps~^. The theory shows qualitative agreement with 
the experiment. The insensitivity of the calculations to 
the frequency dependence of relaxation rate (minor dif- 
ference between solid and dashed lines in Fig. [5]) sug- 
gests that effects of the band structure play the dominant 
role in determining the shape of the infrared conductiv- 
ity and overshadow the strong frequency dependence of 
T obtained within our model and presented in Fig. [5] 

An alternative possible experimental probe that could 
reveal the details of the frequency and momentum de- 
pendence of the carrier relaxation rate in more explicit 
ways are measurements of the position and lineshape of 
the plasmon resonance itself. It was shown in Ref. [2^ 
that these quantities are sensitive to the carrier relax- 
ation time, with both real and imaginary part of r and 
its dynamic nature being essential. Our approach seems 
to fit well to describe such experiments. 

As was mentioned before, the magnetic impurities 
bring localized spins into the system, which interact with 
the itinerant carriers through the p-d exchange interac- 
tion. The fluctuating part of this interaction constitutes 
the spin disorder. The mean field part of exchange inter- 
action, which we absorb into the clean system band struc- 
ture Hamiltonian H^, is responsible for the spin splitting 
of the valence bands once the system enters the magnet- 
ically ordered state. Due to the spin-orbit interaction 
within the valence band, this spin splitting strongly de- 
pends both on the magnitude and direction of the wave 
vector k. 

In Fig. [7] we plot the band structure of ferromagnetic 
Gao.95Mno.05As. Strong anisotropy of the valence band 
spin-splitting is seen between directions along and per- 
pendicular to the magnetization of localized spins {z- 
direction). The inset shows a cut of the Fermi surface 
by the plane ky = 0. One can easily see the distortion of 
the Fermi surface from the spherical shape of the para- 
magnetic system (for clarity we have neglected here the 
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FIG. 9. (Color online) Temperature dependence of the in- 
frared conductivity of Gao.948Mno.052 As with hole concentra- 
tion p = 3 X 10^° cm~^ and Tc — 70 K. Upper panel: experi- 
mental data of Ref. [U- Lower panel: results from Eq. H12p . 



valence band warping, but it is included in our calcula- 
tions). The modification of the Fermi surface together 
with the suppression of localized spin fluctuations are re- 
sponsible for the significant drop in static resistivity of 
GaMnAs during the transition from paramagnetic to fer- 
romagnetic state. This effect was considered beforei ^^i^^ 
Here we point out that the modification of the valence 
band structure during the transition from paramagnetic 
to ferromagnetic state also modifies energies and oscil- 
lator strengths of intervalence band optical transitions 
affecting thus the infrared conductivity as well. To bet- 
ter show the underlying physics of temperature induced 
changes, we plot in Fig. [8] the infrared conductivity for 
the sample parameters of Ref. [s^, but with a small life- 
time broadening of F = 5 meV. In the paramagnetic state 
(solid line) three features can be identified: a strong peak 
around 0.2 eV corresponding to the heavy hole - light hole 
transitions, a smaller peak with a broad shoulder around 
0.4 eV associated with the split off to light hole transi- 
tions and a wide background of split off to heavy hole 
transitions. 

With the temperature going below Tc = 70 K, two 
main phenomena occur. The first is the suppression of 
the high energy shoulder of the split off to light hole 
transitions. The second is the appearance of the tran- 
sitions between the spin-split heavy hole and light hole 
bands and the redistribution of the oscillator strength 
among them. The lowest energy peaks correspond to the 
transitions between spin-split bands. Calculations were 
performed for light linearly polarized in the plane per- 
pendicular to the magnetization. Due to the spin-orbit 
interaction within the valence band, the transitions be- 
tween the spin-split states are optically allowed. The 
additional peak at higher energy corresponds to heavy 
hole-light hole "spin-flip" transitions. As the temper- 
ature goes down, the spin splitting increases and the 



"spin-flip" transitions gain the intensities at the account 
of "spin-conserving" heavy hole-light hole transitions. 

The real GaMnAs samples are much more disordered. 
In Fig. [5] we compare experimental data on infrared con- 
ductivity of Gao.948Mno.o52As from Ref. [13 with calcula- 
tions using our model of Eqs. and PH]) . The large 
disorder induced life-time broadening blankets most of 
the features discussed above. The suppression of the high 
energy shoulder of split off to light hole transitions in the 
ferromagnetic state is seen, however, both on the experi- 
mental and theoretical plots. Overall, for energies above 
the main peak position around 0.2 eV, the calculations 
are in good agreement with the experimental results. 

Note also that unlike in Ref. [ij, our calculations do 
not require incorporation of an impurity band within the 
energy gap to avoid a drop in conductivity around 0.8-1 
eV. At energies below the main peak position the agree- 
ment with the experiment is worse. We should mention, 
however, that this is the region of uit < 1 where our 
calculations are less reliable due to approximate nature 
of the expression ([12]). The self-consistent evaluation of 
Eq. dl]) should be used there instead. Once the frequency 
goes to zero, the static conductivity should more appro- 
priately be calculated using an expression derived from 
the semiclassical Boltzmann equation.— We have inves- 
tigated this regime before^^ to describe the drop in static 
resistivity in the ferromagnetic phase. 



IV. CONCLUSIONS 

We have developed a comprehensive theory of trans- 
port in spin and charge disordered media. The theory is 
based on the equation of motion of the paramagnetic cur- 
rent response function of the disordered system, treats 
disorder and many-body effects on equal footings, and 
combines a k • p based description of the semiconductor 
valence band structure with a full dynamic treatment of 
electron-electron interaction by means of TDDFT. We 
have applied our theory to the specific case of GaMnAs. 

We have shown that the multiband nature and spin- 
orbit interaction within the valence band bring new ef- 
fects for p-doped GaAs as compared to the conventional 
n-type systems. The density and spin-density responses 
of noninteracting carriers within the valence band are not 
the same anymore. Moreover, the long wavelength side 
of the single-particle excitation spectrum is now com- 
pletely dominated by the inter- valence band spin excita- 
tions. Due to the extended region of single-particle ex- 
citations within the valence band, the collective plasmon 
mode entirely falls within the region of these excitations 
and, therefore, is effectively damped for all wave vectors. 

For the magnetically doped system the mean-field part 
of the p-d exchange interaction between itinerant holes 
and localized spins substantially modifies the semicon- 
ductor band structure once the system enters a magnet- 
ically ordered phase. This modification substantially af- 
fects energies and oscillator strengths of the intervalence 
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band optical transitions. Our calculations are in good 
agreement with experimental data for the temperature 
dependence of the infrared conductivity in GaMnAs. 

A full dynamical treatment of electron-electron inter- 
actions is essential to capture the influence of the collec- 
tive excitations on the carrier relaxation rate. Our cal- 
culations show that, by providing an effective channel of 
momentum relaxation, the collective excitations within 
the valence band significantly (up to 50%) increase the 
transport relaxation rate. 

However, it turns out that the actual infrared absorp- 
tion spectra are not very sensitive to the details of the fre- 
quency dependence of the relaxation rate, but are mostly 
determined by the features of the band structure. Direct 
measurements of the position and lineshape of the plas- 
mon resonance itself are likely to be more sensitive to 
the details of the frequency and momentum dependence 
of the carrier relaxation rate. 

The theory presented here, treating disorder and 
many-body effects on equal footings, provides a very gen- 
eral framework for describing electron dynamics in ma- 
terials. It can, in principle, be made self-consistent and 
thus be applied beyond the weak-disorder limit; it can 
accommodate many different types of disorder, as well as 
band structure models. This should make it well suited 
for further exploration of the optical and transport prop- 
erties of DMSs and other systems of practical interest. 



ACKNOWLEDGMENTS 

This work was supported by DOE under Grant No. 
DE-FG02-05ER46213. 



Appendix A: Generalized k ■ p approach 

The derivation of the generalized k • p perturbation 
approach presented here is based on Ref. [23. First, the 
electronic wave function is expanded in the Luttinger- 
Kohn basis'*^ 

* = ^ A„(k)x„,k = ^T. M^)e''"'\un,o), (Al) 

n,k n,k 

where |un,o) are periodic parts of Bloch functions at k = 
0, and A„(k) are the expansion coefhcients. This results 
in the following matrix form of the Schrodinger equation: 



n.k 



Sn.Q H- 



2mo 



— £ 



A, 

mo 



-0, 



(A2) 



where e„_o are the band edge energies at k = 0. 

The last term in Eq. (IA2p mixes states with different 
n for k ^ 0. Now we separate the whole set of the bands 
{n} into those whose contribution we are going to calcu- 
late exactly {s}, and the remote bands {r} that we will 
treat up to the second order in momentum. Equation 



2]) can be represented as 

{Ho + Hi+ i?2)A = eA, 



(A3) 



where A is the vector of coefficients A„(k), Hq is the di- 
agonal part of Hamiltonian, and Hi and H2 correspond 
to the block-diagonal and off-block-diagonal parts of the 
k • TT matrix with respect to the included and remote 
bands. Next, we apply the canonical transformation 



A = e^B = e 



B, 



(A4) 



with 5*1 and 5*2 being antihermitian operators of first and 
second order in the perturbation, respectively. The ma- 
trix equation (jA3[) then has the form 



{Ha + Hi+ H2)e^'+^^}B = HB = eB. 



(A5) 



By choosing 

H2 + [Ho, Si] = 0, [Ho, S2] + [Hi, Si] = 0, (A6) 

where [...] denotes the commutator, we write up to terms 
of second order in the perturbations Hi and H2 



H^Ho + Hi + ^[H2,Si]. 



(A7) 



The matrix elements between the Luttinger-Kohn peri- 
odic amplitudes |u,i.o) = [n) are 



{n[Ho[n') = [SnM + ) Sn,n', 

{s[Hi[s')=Y.^^. 

a mo 

{s[H2[r)^Y.^^^ 
^ Too 



{s[Si[r) 



{s\Ho\r-) 



s[Ho[s) - {r[HoV) 
hk Kr 1 



Too Erfi — £s,0 

For the last term in (IA7I) we can then write 



(A8) 
(A9) 
(AlO) 

(All) 



{s[[H2,Si][s')=J2{{s[H2[r){r[Si[s') 

r 

-{s[Si[r){r[H2[s')] (A12) 



E 



Here we used the fact that the H2 and Si operators have 
only off-block-diagonal matrix elements between the s 
and r bands. Eqs. (jA8P - (jA12[) define the matrix of the ef- 
fective Hamiltonian (IA7|) . Nonvanishing matrix elements 
are determined by the symmetry of the crystal. 
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Appendix B: Evaluation of the matrix elements in Eqs. H17|) and (|18|) 

In order to evaluate Eq. ((T7| we need to calculate the following matrix element: 



where |ui,k) is expressed through the amplitudes at the zone center: 



(Bl) 



(B2) 



From diagonalization of the effective Hamiltonian (jA7|) . however, we obtain coefficients i?„(i, k) related to A„(i, k) 
through Eq. (jA4l) . Expanding e'^ « 1 + 5, we express 



(B3) 



where we have used the fact that the coefficients i?„ are non-zero only for exact bands and S has only off-block-diagonal 
matrix elements. The bra vector is 



(u.^k'| = E^:'(*''k')(svEE(^'|'^(k')kW(^^k')(r'l 



(B4) 



where we have used the antihermiticity of S. Matrix elements of an arbitrary operator F to the lowest order in S can 
then be expressed as follows: 



s's \ r 

Using Eq. (IA11[) for matrix elements of , we have 



M'\F\u,M)^J2^s'i^'X)B,ihk) {s'\F\s) E 



kx{s'\F\r){r\Tr^\s) k'^{s'\7r^\r){r\F\s) 



Too 



(B5) 



(B6) 



The matrix element (|Bip has thus the following form: 



^i' ,k— q 



n (fca - y) + M,,k) = E - q)5s(«, k) 



Too 



S s £r 



For q = it reduces to 

{Ui,,i^\hka +TTa\Ui^u.) = ^ B*, {i'^k) B sii,'k) 



Too 



E^^ 



7r° TT^ n\ tt" 

s' ,r r.s s' .r r.s 



By comparison with the expressions derived in Appendix A, we find that this reduces to 

("i',k \^ka + TTa\Ui,k) = E -^s' (*' > k)^^ (l, k) - 



, mo d /„/|ff|„\ 
Ti dka 



where H is the Hamiltonian (jA7p . 

The matrix elements of the spin operator in Eq. (|18p should also be evaluated through Eq. (|B6 



(B7) 



(B8) 



',k'|<^^Kk> =E^:'(*''k')i?.(*,k) I {s'\a^\s) ~ —Y, 

A.r 



-r "-s' 



(B9) 



Let us look now at the sum over remote bands. Since the spin operator acts only on the spin part of the basis 
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functions, only those remote bands whose orbital part has 
the same symmetry as the exact bands will contribute to 
this sum. 

If we are considering a 6 x 6 Hamiltonian and neglect in- 
version asymmetry, the exact states arc p-bonding states 
that transform according to the representation of 
the point group Oh {T[q small representation). The mo- 
mentum operator transforms as Fj" , and since the direct 
product X F2 X F^ does not contain a unit repre- 
sentation, the sum over remote bands vanishes. There 
may be a small contribution in Td crystals, but it can be 
considered negligible. 

If we are working in an 8-band k • p model, there are 
possible contributions to the sum when \s) and |r) are T[ 
states and |,s') is F'j^g and vice versa. Since there is only 
a small admixture of the conduction band amplitude to 
the valence band states, these contributions are expected 
to be small and therefore can be neglected. 

Based on this reasoning, we use the following approx- 
imation: 

{ui',u'\^''\ui,u) ~ 5^S;,(i',k')B«(i,k)(s'|a''|s)- (BIO) 



Appendix C: 8 x 8 Hamiltonian 

In the basis 



1) = 


1, 


St, 


2) = 


1 

\E,-^) = 
1 > 2' 


iS i, 


0; — 


\HH +-) 


— ^— ( X -i- iV\ t 

- ^(A+zy)T, 


4) = 


\LH,+-) 


= —^[{X + iY) i -2Z tj, 


5) = 


\LH,--) 


= -^[{X - lY) t +2Z i\. 


6) = 


\HH,-l) 




7) = 


\soy-) 


= l=[{X + iY)i+Zt], 


8) = 


\so,-l) 


= ^[-{X-zY)t+Zi], 



the Hamiltonian matrix has the form 



2me 



V2 



Vk+ 
















P + Q 




-^Vk- 


L* 






M* 














* 7- /* 


7f'- 




iV2M'* 



Vk, 



if'* 



P-Q 


M* 
iV2Q' 



^/6 



Vk- 



Ivk. ^Vk. 



M 


P-Q 
-L* 

Is 



2 

iV2Q' 





M 
-L 



Vk, 



\/2 

Pi.,. 



\/3 



-i\f2M' 



2 



iV2Q' 



P + Q -iV2M'* ^L'* 



L' iy/2M' P' -A 



\/2 



P'- A 



(C2) 



with 



k-\- — ^^3/) 

V = -i—{S\p^\X) 
mo 



En 



2mo 



Interaction with remote bands results in the intra valence 
band terms 



pC) = - 

gw = - 



^ ~(');2 



^^\kl + kl-2kl). 



2mo 

-2 



2mo 
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-V3 



2mo 

where renormalization leads to 



1 



71 



1 



1 Er, 



71 = 71 



mo 
Ep 
3^' 



71 



72 = 72 



3(i?g. 

Ep 
12 



72 - TTT TT 



*73 {kxky 



1 

Ea + A 



Er, + A 



73 = 73 



73 = 73 



Ep 

6Eg 

Ep 
T2 



Ea 



This reflects the fact that the interaction between con- 
duction and valence bands is taken in our Hamiltonian 
explicitly. In writing the matrix (jC2p we have neglected 
small terms associated with the lack of inversion symme- 
try in Td crystals. 

The matrix of the mean- field part of the s{p)-d ex- 
change interaction, which is responsible for the band spin 
splitting in the magnetically ordered phase, has the form 



/ 


a 




















\ 







—a 














































































































































V 














i^B 

* V3P 








) 



(C3) 



where the 2;-axis is chosen in the direction of the magne- 
tization and Noa and Nq/S are the s-d and p-d exchange 
constants. 

The mean field value of localized spins is determined 
as the thermodynamical average 



(S) = (S,) = 1 Tre-^5; 

with the partition function 



(C4) 



(C5) 



Within the mean field approximation for uncorrelated 
spins the spin Hamiltonian is 



Hm — —BcsSz, 
with the efi'ective field 

and 

3kTr 



Jo 



S{S + 1)' 



(C6) 
(C7) 

(C8) 



The Curie temperature Tc is an input parameter of our 
model; through the transcendental equations (|C4I) and 
(jC7p it determines the mean field value of (5) . 



Appendix D: Local field factors for partially spin 
polarized systems 

Expressions for local field factors of partially spin po- 
larized electron gas were derived in Ref. '33', but in a 
different spin basis. Here we will briefiy rederive them in 
the basis of Eq. (g]). 

In the adiabatic approximation (which ignores fre- 
quency dependence), the components of the tensor /xc 



(Dl) 



of the local field factors in Eq. ([23| have the form 
52 [nebcin, 0] 



dptOpj 



where Cxc is the xc energy per particle, n = pi is the 
electron density, and ^ is the spin polarization: 



C = 1^ = -\ pI + lip+p- + p-p^ 

n V 2 



(D2) 



We assume here that Cxc depends only on the absolute 
value of |^|. Direct evaluation of Eq. (ID1[) gives 



/ii = 2$^-2e"'^ 



9^exc , ^^Cxc 



dpi ''dpid^ dpi ' pi c?^2 



dp, [^''dpid^ ^ oe 

h.=A + pIB, 
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fz- 


2 


/++ 


- P-P- R 
^ 4 ' 


f— 


" 4 ^' 


f+- 


2 4 



with 



1 ^e^c 



B 



1 9^6x0 Se^c 



Note that /ii' = /^'i and, generally, the tensor of local 
field factors is a symmetric matrix. If, however, the z- 
axis is directed along the average spin direction, so that 
the ground-state transverse spin densities vanish (p+ = 
p_ =0), then the matrix reduces to the block-diagonal 
form of Eq. 

We define the xc energy of the spin polarized system 
in the usual manner as^ 

0)+(e.e(n,l)-e.e(n,0))/(O, (D3) 

with 

n^) 2(21/3-1) • ^^^^ 

This is exact for the exchange part, but only approxi- 
mately so for the correlation part (which will be neglected 



anyway in the following). With this, we get 
i9ex 



(e.e(n, 1) - e.e(n, 0)) ^ 3(2i/3 1 i) ' 



(^exc(n, l)-exc(?^,0) 



(21/3 _ 1) 

(D5) 

(l+g)-2/3^(l_g)-2/3 
1(21/3 - 1) 



and 

(D6) 

This completes the definition of the local field factors 
for partially spin polarized system. The only remaining 
ingredients we need to perform the actual calculations 
are the expressions for the xc energy for unpolarized and 
fully spin polarized system, exc('T-, 0) and exc(?^, !)■ In 
this work for simplicity we limit ourselves to the exchange 
part of Cxc: 



ex(n,0) = -— 1^- 
ex(n,l) = 2i/3ex(n,0), 



(D7) 
(D8) 



where K is the static dielectric constant of the host ma- 
terial. Direct evaluation gives the following expressions: 

1/3 



9exc 


e2 


dn 


8K 


5^exc 


e2 


dn? 


12K 


5^exc 


e2 


dn 9^ 


6K 



5/3 ((1 + ^)4/3 ^(1^^)4/3^^ 



1/3 
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